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Final  Report 

VISCOUS  EFFECT  ON  SHIP  WAVE  RESISTANCE 

Work  on  the  subject  problem  under  the  sponsorship  of  DWTNSRDC 
was  initiated  in  October  1981  and  terminated  in  September  1982.  This 
research  is  being  continued  under  the  Special  Focus  Program  of  the  Office 
of  Naval  Research* 

During  the  contract  year,  one  paper  was  published  and  two  reports 
were  produced.  These  are  the  following: 

1.  "Centerline  Distribution  as  Solutions  of  Neumann  Problems  for 

Ogivelike  Forms”f  by  L.  Land weber,  Eighteenth  Midwestern 

Mechanics  Conference,  May  1983. 

2.  Progress  Reports  by  Kazuhiro  Mori. 

A.  Calculation  of  Wave  Resistance  and  Sinkage  by  Rankine- 
Source  Method. 

B.  Prediction  of  2-D  Near  Wake  Flow  by  Making  Use  of  Time- 
Dependent  Vorticity  Transport  Equation. 

C.  Free-Surface  Boundary  Layer  and  Necklace  Vortex 
Formation. 

EHR  Report  No.  262,  The  University  of  Iowa,  May  1983. 


3. 


"Numerical  Evaluation  of  the  Havelock  Green  Function”,  by 
Lawrence  K.  Forbes,  unpublished.  This  is  attached  as  part  of  this 


CENTERLINE  DISTRIBUTIONS  AS  SOLUTIONS  OF  NEUMANN  PROBLEMS 

FOR  OGIVELIKE  FORMS 


L.  Landweber,  Iowa  Institute  of  Hydraulic  Research 
The  University  of  Iowa,  Iowa  City,  Iowa  52242 


Introduction 

As  in  well  known,  approximate  solutions  of  symmetric  Neumann  problems  for 
thin,  symmetric  two-dimensional  forms  can  be  obtained  by  formulating  the 
boundary  condition  as  a  Fredholm  integral  equation  of  the  first  kind.  Since 
exact  solutions  of  such  integral  equations  exist  only  when  certain  conditions 
are  satisfied.  It  Is  of  Interest  to  find  cases  for  which  exact  solutions  can 
be  obtained. 

The  current  Interest  in  this  problem  stems  from  an  attempt  to  find  a 
centerplane  distribution  for  a  slender  ship  form  with  ogive-like  sections. 
Slender-body  theory  would  yield  an  approximate  centerplane  distribution  for 
such  a  form  if  the  particular  Neumann  problems  posed  by  the  slender-body 
theory  have  exact  solutions  for  the  transverse  sections  of  the  double  ship 
form. 

In  the  present  work,  we  shall  be  mainly  concerned  with  solutions  for 
ogival  sections.  Applications  of  the  procedures  to  other  forms,  such  as 
elliptical  and  doubly-parabollc  (Wlgley)  sections  will  be  indicated. 

Geometry 

The  ogive  is  a  two-dimensional  form  derived  from  the  intersection  of  two 
circular  arcs  of  equal  curvature.  For  a  form  of  thickness  2c  and  length  2h, 
its  equation  Is 

M  *  7c  +  c2)2  "  4c2j<2]  ^  -  h2  +  c2} ,  -  h  <  x  <  h  (1) 


(2) 


its  slope  at  x  *  -  h  is 

tan  a  -  -J—%  .  c  =  h  tan  | 

and  Its  radius  of  curvature  is  R  *  h  sec  a.  The  cosine  of  the  angle  of  the 
normal  to  the  ogive  with  the  x-  axis  is 


3x 

3n 
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(3) 


The  exterior  of  the  ogive  in  the  z  =  x  +  iy  plane  is  mapped  into  the 
exterior  of  the  unit  circle  in  the  plane  c  =  C  +  i  n  *  p  e^  by  the  successive 
transformations 


where 


T 


-  1 


(4) 


1  x  <  1  (5) 

Mappings  of  the  ogive  in  the  z,  Z,  T  and  planes  are  shown  in  Fig.  1.  The 
upper  half  of  the  ogive  maps  into  the  crescent  ACBDA  in  the  c-  plane;  the 
lower  half  (not  shown)  would  give  the  mirror  image  of  the  crescent. 

Let  e2,  <|> ,  $2  denote  the  polar  angles  in  the  Z,  5,  and  T-  planes,  and  put 


u  *  (tan  $/2)  ,  E  *  [  m2  +  2o>  cos  a  +  1  ]  ^2  (6) 

Then  $2  3  *  ®2  and  the  ^ocat^ons  points  in  the  mappings  of  Fig.  1  are  given 
in  the  following  table: 
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Solutions  of  the  Neumann  Problem 


We  wish  to  find  solutions  of  the  Laplace  equation 


32« 

7? 


(7) 


for  the  exterior  of  the  ogive  which  have  zero  circulation 
satisfy  the  boundary  condition 


3« 

SIT 


f(x) 


and,  on  its  contour. 


(8) 


where  n  denotes  distance  normal  to  the  ogive,  positive  outwards.  The  form  of 
the  Neumann  boundary  condition  (8)  implies  that  3$/3 n  is  the  same  at  ±  y(x). 
The  flux  through  the  ogive  is  given  by 


Q  =  2nM  =  2  J  f(x)  dx  (9) 

-  h 

where  s  denotes  arc  length  along  the  ogive,  positive  in  the  counterclockwise 


sense. 


In  the  t-  plane,  the  Neumann  boundary  condition  becomes 


s  f<x>  1^1  ’  o  a  1 


For  the  ogive,  we  have 


x  =  h  (1  -  a)2)/  E2 


I  del - 12  - 

^  X  E  sin  4> 

We  define  a  zero-flux  potential  by 

♦  -♦.♦H  (13) 

The  complex  potential  for  this  Neumann  problem  is  then 

2tt 

w(c )  *  *  +  i*  =  M  In  c  +  ^  /  |~r  In  (c  -  e1*  )  d(j> ' ,  p  >  1  (14) 

and  the  complex  velocity  Is 

dw  -  M  .  1  2*  3«  <%'  s  ,  nn 

The  form  of  dw/d?  in  (15)  cannot  be  used  to  continue  the  complex  velocity 
analytically  into  the  unit  circle.  If  none  of  the  singularities  of  dw/d?  lies 
within  the  crescent  ACBOA  in  Fig.  1,  the  function  can  be  continued 

analytically  to  the  arc  BCA  In  the  plane,  and  hence  to  the  axis  of  the 

ogive  in  the  z-  plane.  We  would  then  have,  at  the  ogive  axis, 


a  ff  %  "  u  (x’°)  -  iv  (x*  0+) 


(16) 


where  u  and  v  are  velocity  components  in  the  x-  and  y-  directions,  and  0+ 
indicates  that  y  ♦  0  through  positive  values.  The  axial  source  distribution 
would  then  be  given  by 


m(x)!-v  (x.O)  *  -  “  I" 


where  Im  denotes  the  imaginary  part. 

Exanple:  Ogive  In  a  Uniform  Stream 

For  this  case,  the  exact  centerline  distribution  can  be  found  by  a 
simpler  method  and  is  given  by 


■(*)  .  JJ  xVsIn  2.X  V*  >  -  — 4T7  (W) 

[(h  +  x),x  -  2 (hz  -  xzr*  cos  2.X  +  (h  -  x)4^ 

More  generally,  however,  an  explicit  algebraic  formula  for  the  complex 
velocity,  as  was  used  to  derive  (18),  is  not  available.  An  alternative 

procedure  is  to  investigate  the  singularities  of  the  integrand  of  (15), 

expressed  as  a  line  integral  around  the  unit  circle  of  a  function  of 

C 1  -  e  .  If  3$/3p'  can  be  expressed  as  a  function  F(c'),  (15)  could  be 

written  as 

37 * ? * F(c') n  (19) 


.4X....2  2.  2x  -  1 

h  -  x)l(h  -  x  ) _ 


In  the  present  case,  M  =  0  and,  by  (3) 

,  3*  3 x  x 

'(*>  =  5F  ■  -  sr  *  -  K 

and,  by  (10),  (11)  and  (12), 


sm  a 


3#  _  2h(u  (1  -  hi  )  sin  a 


4“ 

X  E  sin  $ 


‘  ,  p  '  =  1 


Since  sin  $'*(?'  -  l)/2i?  ,  and  one  can  show  that 


-  -  -  (rrv 
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we  see  that  3*/3p'  can  be  expressed  as  F(c'),  and  (19)  becomes 


dw  .  4  h  sin  q 


i,  hi  (1  -  hi  )  d$  1  v  , 

P  -T -  p - - -  ,  P  >  1 

E4  fc'*  -  1)(C  -  C  * ) 


The  apparent  singularity  at  c ’  *  ±  1  is  removable.  The  partial-fraction  form 


]  (24) 


M  (1  -  M2)  ,  1  r _ 1 _ 1  _ 1_  +  1 

V  2  sin  a  Lwe-io+1  {ue-i<*+1)2  '  ueia+l  (ajei«+1 )2 

then  shows  where  poles  of  the  integrand  may  be  present.  We  find  that 
(we'io+l)  has  no  zeros  within  the  unit  circle,  but  that 


wtia  +  1  =  2  (pc'  -  uV2  +  ...)  (25) 

vanishes  at?'  *  0.  Hence,  since  there  are  no  singularities  within  the 
crescent,  this  confirms  the  existence  of  an  exact  centerline  distribution. 

It  is  also  of  interest  to  attempt  to  evaluate  the  integral  in  (23).  In 
order  to  apply  the  Cauchy  residue  theorem,  it  is  necessary  that  the  integrand 
be  a  regular  function,  except  at  g '  =0.  Because  of  symmetry,  the  imaginary 
parts  of  the  integrands  at  conjugate  points  in  the  unit  circle  must  have 
opposite  signs  when  g  *  £  >  1.  These  two  requirements  can  be  satisfied  only  if 
the  integrand  is  real  along  the  diameter  g'*C',-l<5'  <1.  It  can  be 
shown  that  this  last  condition  is  satisfied  by  the  third  and  fourth  terms  of 
(24).  Since  these  are  also  the  terms  which,  by  (25),  yield  the  singularity, 
one  may  apply  the  residue  theorem  to  that  part  of  the  integrand  containing 
these  two  terms.  By  applying  (25),  we  obtain 

_  i  r  i  l  i  _  l  , 

U  *  C  )\5  -  1)  (we  +1)  we  +1  4  i j  c  C 

The  expression  (23)  for  the  complex  velocity  then  becomes 


dw  _  2ih  r  r  1 _ 1  _ dg 1  _  Xjh 

*  "X  9  we'  ia+l  "  (we"ia+lr  (c  ,2-l)(c-c  ’ )  ’  C2 

Comparison  with  the  known  exact  expression 


dw 


»» (i  •  V  •  I 

c 


(26) 


(27) 


indicates  the  value  of  the  remaining  integral  in  (26).  Since  the  present 
example  was  treated  for  the  purpose  of  illustrating  a  procedure  for  finding 
centerline  distributions,  and,  in  general,  integrals  such  as  that  in  (26), 
would  require  numerical  evaluation,  further  consideration  of  this  case  would 
not  be  rewarding. 
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Procedure  for  Ogivellke  Sections 

1.  Map  the  given  section  in  the  z-  plane  into  the  unit  cirt/cle  in 
the  c-  plane.  For  an  ogivelike  section,  this  is  best  done  with  the 
preliminary  transformation  (4)  to  transform  the  section  first  into  a  nearly 
circular  form. 

2.  Map  the  axis  of  the  ogive  into  the  interior  of  the  unit  circle.  The 
crescentlike  region  bounded  by  this  curve  and  the  arc  of  the  unit  circle  gives 
the  mapping  of  the  upper  (or  lower)  half  of  the  given  section. 

3.  Obtain  3*/3p,  the  Neumann  boundary  condition  on  the  unit  circle, 
express  it  as  function  F{c),  c  =  e1*,  and  determine  the  singularities  of  the 
integral  expression  (15)  for  dw/dc.  If  none  of  the  singularities  lies  within 
the  crescentlike  region,  then  there  exists  a  centerline  distribution  of 
sources  which  gives  the  solution  of  the  given  Neumann  problem  for  the  section. 

4.  If  the  solution  exists,  it  can  be  obtained  by  solving  numerically  a 
Fredholm  integral  equation  of  the  first  kind  formulated  by  equating  the  normal 
component  at  the  surface  of  the  given  profile,  induced  by  the  unknown  axial 
source  distribution,  to  the  prescribed  Neumann  boundary  condition. 
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NUMERICAL  EVALUATION  OF  THE  HAVELOCK  GREEN  FUNCTION 


by 

Lawrence  K.  Forbes 


Iowa  Institute  of  Hydraulic  Research 
The  University  of  Iowa 
Iowa  City,  Iowa  52242  USA 


Abstract 

Three  closed-form  approximations  to  the  double-integral  "near-field"  term 
in  the  Havelock  Green  function  are  derived.  An  application  to  centerplane 
source  distributions  for  ship-wave  problems  is  discussed. 
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1.  Introduction 

In  the  theory  of  potential  flow  about  arbitrary  bodies  moving  steadily  in 
or  beneath  a  linearized  free  surface,  the  velocity  potential  is  usually  sought 
as  the  solution  to  an  integral  equation  having  as  the  kernel  a  derivative  of 
the  Havelock  Green  function  in  the  direction  normal  to  the  body  surface.  The 
choice  of  the  Havelock  Green  function  as  the  fundamental  solution  to  Laplace's 
equation  is  natural,  since  it  immediately  ensures  that  the  linearized  free- 
surface  condition  and  the  radiation  condition  are  satisfied.  However,  this 
important  advantage  is  offset  by  the  great  complexity  of  the  function  itself, 
forcing  some  investigators  [1,  2,  3]  to  abandon  it  completely  in  favor  of  the 
simpler  Rankine  source  function. 

Because  of  the  central  role  played  by  the  Havelock  function  in  the  theory 
of  ship  wave  resistance,  its  numerical  evaluation  has  been  the  subject  of  much 
research.  There  are  many  different  expressions  for  this  function,  all  of 
which  may  be  related  to  one  of  the  three  basic  forms  summarized  by  Noblesse 
[4],  who  simplified  the  form  of  the  function  somewhat  through  the  introduction 
of  the  complex  exponential  integral. 

In  the  numerical  application  of  the  Havelock  Green  function  to  the 
solution  of  practical  ship-wave  problems,  there  are  two  basic  difficulties  to 
overcome.  The  first  of  these  is  due  to  the  additional  singular  behavior  of 
this  Green  function  on  the  plane  of  the  linearized  free  surface,  described  by 
Noblesse  [4],  This  problem  may  be  circumvented  by  adopting  the  approach  of 
Miloh  and  Landweber  [5],  They  showed  that,  under  certain  assumptions  about 
the  ship-hull  geometry,  the  problem  could  be  reduced  to  finding  an  appropriate 
source  distribution  M  on  the  ship  centerplane  satisfying  the  f i rst-kind 
Fredholm  integral  equation 


...  Jl.. 


3 


r. 

E; 


s 

r? 

I 


si 


1 


/  M  (P)  dSp  +  §jj-.  0  (1) 

on  the  ship  hull 

y  =  ±  f ( x , z)  ,  (2) 

where  P (c »0 )  and  Q(x,y,z)  are  points  on  the  centerplane  and  ship  hull, 
respectively.  The  Havelock  Green  function  is  denoted  by  the  symbol  G(Q,P)  and 
Nq  represents  the  direction  normal  to  the  hull.  By  discretizing  the  integral 
on  the  left-hand  side  of  equation  (1)  using  Gaussian  quadrature,  the 
requirement  that  numerical  grid  points  be  placed  on  the  plane  z  =  0  of  the 
undisturbed  surface  is  eliminated,  and  the  singular  nature  of  the  Havelock 
Green  function  on  this  plane  is  avoided.  The  results  of  the  present  paper 
have  been  developed  primarily  with  a  view  to  their  eventual  implementation  in 
the  numerical  solution  of  equation  (1). 

The  other  major  difficulty  with  the  Havelock  Green  function  is  the 
computing  time  required  to  evaluate  it.  On  the  ship  hull,  in  particular,  the 
"near-field"  component  of  this  Green  function,  which  may  be  expressed  in  terms 
of  a  double  integral,  cannot  reasonably  be  ignored  and  an  efficient  method  for 
its  evaluation  must  be  developed.  A  simple  calculation  serves  to  illustrate 
the  point;  if  a  solution  to  equation  (1)  is  sought,  employing  50  points  from 
bow  to  stern  and  20  points  from  keel  to  free  surface  in  the  discretization  of 
the  integral  (a  total  of  1000  points),  then  10®  separate  evaluations  of  the 
Green  function  would  be  required  to  form  the  kernel  of  the  integral 
equation.  A  single  accurate  evaluation  typically  requires  in  excess  of 
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1/10  th  sec.  computing  time,  so  that  approximately  10  secs,  or  one  day 


continuous  computing  time  would  be  required  to  solve  the  problem,  which  is 
unreasonable.  Such  considerations  lead  Lee  [6]  to  conclude  that  the  function 
cannot  be  applied  to  the  solution  of  practical  flow  problems. 

In  the  present  paper,  a  number  of  different  approximations  to  the  basic 
function  ([7],  page  484) 
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and  its  derivative  in  the  direction  normal  to  the  hull  (2)  will  be  derived, 
and  their  suitability  for  the  rapid  evaluation  of  these  functions 
determined.  Here,  F  is  a  Froude  number  based  on  the  ship  length  and  R  is  the 
distance 

R  =  [x'2  +  y'2  +  (z  -  C)2]  ^2  (4) 

between  the  points  Q(x,  y,  z)  and  P(£,  n,  c).  The  quantities  x',  y 1 ,  z'  are 
defined  by  the  relations 

x'  =  x  -  5 

y'  ■  y  -  n  (5) 

z'  =  z  +  ?. 


In  addition. 
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R' 


+  z 


(6) 


and 


p  =  [x*2  +  y 1 2  3 1/2  (7) 

a  =  arctan  (y'/x1 ) 

2.  Direct  Numerical  Evaluation 

In  order  to  evaluate  the  Havelock  Green  function  numerically,  the  Cauchy 
principal-valued  integral  appearing  in  equation  (3)  is  first  transformed  by  a 
straightforward  contour  integration  due  to  Havelock  ([8],  page  291).  After 
some  algebra,  equation  (3)  becomes 
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where  the  functions  R,  R',  etc.  are  defined  by  equations  (4)  -  (7). 

The  (indefinite)  integral  with  respect  to  k  in  equation  (8)  may  be 
evaluated  using  the  Laguerre  polynomial  approach  (see  [9],  page  890, 
formula  25»4«45).  To  demonstrate  that  acceptable  accuracy  could  be  maintained 
by  the  Laguerre  polynomial  method  when  the  integrand  is  oscillatory,  the  test 
functions 


6 


J 

0 


-  x 


,COS  X, 

^in  x* 


dx 


were  evaluated  using  the  fifteenth-order  integration  formula,  and  it  was  found 
that  the  result  agreed  with  the  exact  answer  1/2  to  six  figures  in  both  cases. 

Once  the  k-integration  in  the  double-integral  term  in  equation  (8)  has 
been  performed,  it  remains  to  integrate  the  result  with  respect  toe.  This 
may  be  achieved  using  Gaussian  quadrature.  Similarly,  the  single-integral 
term  in  equation  (8)  may  be  evaluated  by  Gaussian  integration. 

Finally,  the  normal  derivative 


as  {Ml 

9Nq 


»  VpG(Q,P)  -  N 


(9) 


is  usually  required,  as  in  equation  (1),  for  example.  Here, 


N  = 


’  fx  1  +  '  fzk 

a  ♦  t\  +  f zJ/z 


(10) 


is  the  normal  vector  to  the  ship  hull.  The  derivatives  of  G  in  the  x  -,  y  - 
and  z-directions  are  computed  using  forward  differences. 

3.  Asymptotic  Expansion  for  Small  Froude  Number 

To  develop  an  asymptotic  expansion  for  small  Froude  number,  equation  (3) 
is  first  expressed  in  the  form 


G  =  -  ^  i  +  ^3  +  ^3 


y 


where 


.**  »’*  «.**  «*'  •*•>*•*  •  '*  «'*  _%  ,.**  .*•  , 


-v.  ^ 


and 


TT  +  ^7  Re 


Ex  (Z3)  dt 


U3  =  -  -4  Im  ;’/2  sec2e  exp  [*'  8  *  ^  ‘f*  9  cos  <e  -°>]  de  (11) 

F  -  ir/2  +  a  F 

The  complex  number  Z3  in  equation  (11)  is  defined  as 
Z3  *  J2  [z<  (1  +  t2)  +  i  (  x*  +  ty')(l  +  t2)1^  , 
and  Ej  (Z3)  is  the  complex  exponential  integral 

E1  <Z3>  *  /"  • 

z3 

Equation  (11)  is  the  third  of  the  forms  of  the  Havelock  Green  function  derived 
by  Noblesse  [4], 

Since  the  evaluation  of  the  exponential  integral  of  the  complex  argument 
Z3  is  apparently  more  time  consuming  than  the  corresponding  Laguerre- 
polynomial  integration  in  section  2,  equation  (11)  does  not  appear  to  be  of 
much  value  in  the  direct  numerical  evaluation  of  the  Havelock  Green  function, 
despite  the  fact  that  the  double-integral  term  has  been  reduced  to  a  single 
integral  involving  the  exponential  integral,  which  in  some  sense  may  be 
regarded  as  an  "elementary"  function.  However,  equation  (11)  readily  yields 
an  asymptotic  expansion  for  low  Froude  number. 

Since  | Z3 j  becomes  large  as  F  0,  provided  that  x',  y',  and  z'  do  not 
all  simultaneously  vanish,  the  function  £1(^3)  may  be  replaced  by  its 
asymptotic  expansion  for  large  | Z3 | , 


T  ~  »  :■ 


r?1;  T-.  r 


8 


h  (z3)  -  e 


z 

n=l 


(-.Un.  n 

7n+l 

o 


Substituting  this  result  into  the  formula  for  the  near-field  disturbance  term 
N3  in  equation  (11)  yields  the  asymptotic  expansion 


N3~F- 


nr 


Re[/  £ 


z 

n=l 


(-1)  n!  / 


dt 


,n+l 


(12) 


Each  of  the  integrals  in  equation  (12)  may  be  evaluated  in  closed  form,  using 
the  calculus  of  residues,  but  with  greatly  increased  labor  for  every  increase 
in  the  value  of  n. 

We  have  computed  the  first  three  terms  in  the  asymptotic  expansion  (12) 
for  N3,  and  present  the  results  below.  The  first  integral  cancels  with  the 
first  term  on  the  right  hand  side  of  equation  (12);  then,  after  substantial 
calculation,  the  desired  expansion  is  obtained  in  the  form 


2  rx'2  -  y'2  7'  (x'2  -  y'2)  x'2z't 

n3  -  2 r  [x  4  *  +  7  4  +  *3-?-2-] 

p  R  p  R  p 

c  4  rx‘8z‘6  ,  x'6  z'4  (5y ,2  +  z'2) 

•6F  — PT7 

.  »'4  z'2  Illl'V2  *  2z ' 4  »  15v'4)  .  2z‘  (x|4  -  6«i2v'2  +  »'4) 

R'rV  p8 

2  2  4  4 

R'  x,2y'2  (2 1 y '  z'  +  12z'  +  5y'  )  .  R'3  y,4(y'2+  2z'2)  t 

4  8  +  4  8  J 

r  p  r  p 


+  0(F  )  • 


(13) 


where 


To  derive  an  asymptotic  expansion  for  the  wave- like  contribution  W3  to 

the  Green  function  in  equation  (11)  we  note  that,  for  p  *  0,  the 
2 

parameter  p/F  becomes  large  as  F  ♦  0,  and  consequently  W3  may  be  approximated 
using  the  method  of  stationary  phase.  The  result  is  classical,  and  appears  in 
numerous  texts  and  articles  on  ship  waves  (see,  for  example.  Stoker  [10]  and 
Newman  [11]).  This  gives 


^3*  ’  F  .Z-  ^ 


1*1  p  1 9  "  (e  i ) 


u  2  z'  sec  p 9 (0 -; )  pg(ei) 

)/2sec  exp  [ - ^ - ”](cos  ( - 2 — )+sin  ( — 2 — ^  (■ 


p  9  (e  i ) 

V 

F‘ 


pg(ej 

/  _  1  ' 

* 

F< 


-3/2 


for  tan  a  <  2  ,  where  9^  and  9 ^  are 


found  from 


tan  e 


1.2 


:,1. Jl  (1.  -  .8  tan  ff), 


2  te 


4  tan  a 


2 

and  g  (e^ )  =  sec  ei  cos  (e ^ -  a),  i  *  1,  2  . 
When  |tan  a|  >  2  ~^2, 


W3  -  0  ,  (14b) 

-3  /2 

and  for  tan  a  =  +  2  , 

U3  -  -2  3/2  r(y)  ( — 1)1/3  exp  [^]  sin  P^].  (14c) 

pF*  2F^  2F^ 

The  low  Froude  number  form  of  the  Havelock  Green  function  is  thus  given 


by  equations  (11),  (13)  and  (14).  The  normal  derivative  (10)  is  computed 
using  forward  differences. 


9 


roximate  Green  function  for  centerplane  source  distribution  on  a  thin 


An  alternative  form  of  the  Havelock  Green  function  (3),  which  also 
involves  the  exponential  integral  function  Ej,  is  given  by  the  expression 


G  =  "  £  +  ir +  N1  +  W1  • 


where 


2  i  1  Z. 

ij  =  j2  [  »  /  l  lm  (e  Ej  (Zx)  +  In  Zj  +  y}  dt  -  (1 


Wj  =  4H(x')  /  exp  V 


sin  [(-  ty'  -  4)(!  +  dt.  (15) 


The  complex  number  is  defined  as 

Zj  =  t1.  Ci  | x ’  |  -  ty'  +  (1  -t2)1/2z’]  , 


and  H(x')  is  the  Heaviside  unit  step  function,  having  the  value  0  when  x'<  0 
and  1  when  x'>  0.  The  term  In  1\  +  y,  where  y  is  Euler's  constant,  has  been 
included  to  ensure  that  the  integrand  of  the  expression  for  Nj  is  well  behaved 
when  |Zj |  ♦  0. 

Equation  (15)  is  the  first  of  the  forms  of  the  Havelock  Green  function 
derived  by  Noblesse  [4],  who  recommends  its  use  in  the  direct  numerical 
evaluation  of  this  function  (see  also  Noblesse  [12,  13]).  Actually,  equation 


11 


(8)  appears  to  be  better  suited  to  this  purpose,  for  the  reason  outlined  in 
section  3,  In  addition,  the  use  of  equation  (15)  in  an  integral  equation, 
such  as  equation  (1),  requires  special  treatment  of  the  step  discontinuity  at 
x 1  =  0. 

Equation  (15)  may  be  used  to  derive  an  approximation  to  the  normal 
derivative  of  the  Green  function  for  centerplane  source  distributions.  In 
this  case,  the  point  P  lies  on  the  centerplane  n  =  0,  and  the  point  Q  is  on 
the  hull  y  *  f(x,  z).  The  flow  is  assumed  to  be  symmetrical  about  the  plane 
y  =  0,  so  that  only  one  of  the  hull  surfaces  need  be  considered  in  equation 
(2).  If  3 N ^ /a N  denotes  the  normal  derivative  of  the  near-field  term  Nj  in 
equation  (15),  then 


The  argument  of  the  exponential  integral  in  equation  (16)  generally  takes 
moderate  values,  and  accordingly,  the  exponential  integral  may  be  represented 
by  the  formula 


-  zi  e  18  e 

E1  <Z1>  "  e  Cr  +  2  1+6“^  (17) 

1  1  Lx  .x2  +  6. 

This  approximation  Is  due  to  Hershey  [14],  and  clearly  consists  of 
approximating  the  exponential  integral  by  a  series  of  poles  along  the  negative 


real  axis.  Numerical  values  for  e.,  j  =  1,  .  18  and  6  •,  j  *  2,  ....  18 

J  J 

are  given  in  [14], 


Using  equation  (17),  the  first  term  on  the  right-hand  side  of  equation 


(16)  may  be  expressed  in  the  form 


e,  3ir/2  T.f  sgn  x'  -  x'  T? 

-x —  J  a - x - cos  0  de 


irF  -  ir/2 


Tj  +  x' 


■-7I  J 

j*2  J 


where 


3ir/2  (T.  cos  0  +  F  5  .)  f  sgn  x'  -  x'  T?  cos  0  2 

— = - ^ x - x - *-= - cos  0  d0 


(T.  cos  0  +  F  6  .)  +  x*  cos  0 
•  J 


f  sin  0  +  z'  cos  0 


T2  3  sin  0  +  f z  cos  0  . 


The  first  term  in  (18)  may  be  evaluated  immediately,  using  the  calculus 
of  residues,  and  becomes 

2e1  f  z'  sgn  x'  z'  f  +  f  (x'2  +  f2  +  |x'|  R') 


.  .  R'  (|x'|  +  R*) 

The  second  term  could  likewise  be  evaluated  by  the  calculus  of  residues,  by 
summing  appropriate  contributions  from  each  of  the  poles  of  the  integrand,  at 
which  the  denominator  vanishes.  In  this  case,  however,  the  denominator  is  a 
quartic,  and  the  formulas  for  its  zeros  are  too  complicated  for  practical  use. 

If  the  ship  is  very  thin,  f  Is  a  small  quantity  and  terms  involving  it  in 
the  denominator  of  the  second  Integral  in  (18)  may  be  neglected.  Consequently 
the  denominator  becomes  approximately 

z'  cos*  0  +  (2F*6 ,z'  +  x' j  cos  0  +  F  6  -  , 


v 


v 


r 


14 


and 

B  *  (x*2  +  4F26 ,z')1/2.  (19d) 

J 

If  the  argument  of  the  radical  in  equation  (19d)  is  negative,  B  becomes 
imaginary  and  U  and  V  in  equation  (19c)  become  complex  conjugates.  In  this 
case,  Ij  and  l£  in  equation  (19b)  remain  real,  and  (19a)  continues  to  yield  a 
real  result.  In  addition,  a  meaningful  limiting  result  may  be  derived 
for  B  0. 

The  normal  derivative  of  the  wave- 1  ike  term  Wj  in  equation  (15)  is 
obtained  by  straightforward  differentiation.  The  resulting  integral  may  be 
evaluated  by  the  Hermite  polynomial  method  (see  [9],  page  890, 
formula  25»4«46). 


-  ^  /  exp  [l1-  '■  g 
F  0  F 


t  (1  -  t'V*  |y  1 1  r 


-]  cos  [- 


z2  *  ~7  Cl|y'l  l1  +  t2)1/2+  x'  -  tz']. 


This  form  Is  clearly  of  little  use  In  the  direct  evaluation  of  the  Havelock 
Green  function,  since  the  remarks  of  section  3  apply  here  also.  In  addition, 
the  wave- like  part  of  the  Green  function,  W2,  now  contains  two  terms  instead 


of  one. 


An  approximate  form  of  the  Havelock  Green  function,  based  on  equation 
(20),  may  be  derived  as  in  section  4  by  substituting  the  approximate  form  (17) 
of  the  exponential  integral  into  the  expression  for  Ng.  This  yields 


1  +  tfc)  +  (x 


"  j-2 


y‘  t  (1  +  f)  +  ( x 1 1  -  z't  +  F “fij)' 


The  first  integral  in  this  expression  may  be  evaluated  easily  by  the 
calculus  of  residues  and  becomes  simply  ir/(|y'|R').  The  second  integral 
could  likewise  be  evaluated  using  residue  theory,  but  this  would  involve  the 
determination  of  the  roots  of  the  denominator,  which  is  a  quartic  in  t.  This 
is  too  complicated  for  practical  use.  Instead,  the  term  t  (1  +  t  )  in  the 
denominator  of  the  second  integral  in  equation  (21)  is  replaced  by  t  ,  which 
is  exact  for  t  s  0  and  accurate  for  large  t.  Then 


16 


2el  2  18 

-nT  +  ;  I 

R  11  j«2 


i  y1  t2  +  x't  -  z't2  +  6  . 

<1 


The  integral  in  this  expression  may  now  easily  be  evaluated  using  residue 
theory,  and  thus  N£  is  given  approximately  by 


with 


Tj  =  x'2  *  4F2  JjZ' 
t2  =  4F2  4j  |y '  |  • 


(22) 


The  wave-making  term  W2  in  equation  (20)  may  be  evaluated  using  the 
Hermite  polynomial  method  to  perform  the  first  integration,  as  in  section  4. 
The  second  integral  is  evaluated  using  a  Gaussian  quadrature  formula. 


6.  Conclusion 

The  application  of  the  Havelock  Green  function  to  linearized  ship- 

hydrodynamic  problems  typically  results  in  an  integral  equation  to  be 
satisfied  on  the  ship  hull,  with  the  normal  derivative  of  the  Havelock 
function  as  the  kernel.  In  this  case,  the  "near  field"  double-integral  term 
may  be  expected  to  be  of  the  same  order  of  magnitude  as  the  single-integral 
term  responsible  for  the  generation  of  waves  far  downstream,  and  cannot 

meaningfully  be  ignored  in  the  evaluation  of  the  kernel.  However,  the  direct 

numerical  evaluation  of  this  double-integral  term  for  use  in  integral 


equations  appears  to  require  an  unreasonably  large  amount  of  time  on  most 
present-day  computers  [6], 


In  order  to  avoid  the  direct  numerical  evaluation  of  this  douole-integral 
term,  a  number  of  closed-form  approximations  to  it  have  been  developed  in  this 
paper.  Work  is  currently  in  progress  to  utilize  them  to  determine  a 
centerplane  source  distribution  valid  for  the  description  of  flow  about  a 
Wigley  hull  form  at  nonzero  Froude  number,  by  an  iteration  process  similar  to 
that  used  by  Mil  oh  and  Landweber  [5]. 
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Kings  Point,  L.I.,  N.Y.  11024 

The  Pennsylvania  State  University 
Applied  Research  Laboratory  (Library) 
P.  0.  Box  30 
State  College,  PA  16801 

Dr.  B.  Parkin,  Director 
Garfield  Thomas  Water  Tunnel 
Applied  Research  Laboratory 
P.0.  Box  30 
State  College,  Pa 

Bolt,  Beranek  &  Newman  (Library) 

50  Moulton  Street 
Cambridge,  MA  02138 

Bethlehem  Steel  Corporation 
25  Broadway 

New  York,  New  York  10004 
Attn:  Library  -  Shipbuilding 

Cambridge  Acoustical  Associates,  Inc. 
54  Rindqe  Ave  Extension 
Camb r i dqe ,  MA  02140 
\ 

R  &  D  Manager 
Electric  Boat  Division 
General  Dynamics  Corporation 
Groton,  Conn  06340 

Gibbs  4  Cox,  Inc. (Tech  Info  Control) 

21  West  Street 

New  York,  New  York  10006 


Mare  Island  Nav^rfbhipyard 
Shipyard  TecKXcal  Library(202.3) 
Vallejo,  cy'94^98^ 
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Hydronautics,  Inc. (Library ) 
P indell  School  Rd 
Laurel,  Md  20810 


Newport  News  Shipbuilding  and 
Dry  Dock  Company (Tech  Library) 
4101  Wash i ng+on  Aven  ue 
Newport  News,  VA  23607 


! 


/ 


Mr.  S.  Spangler 

Nielsen  Engineering  &  Researc,  Inc. 
510  Clyde  Avenue 
Mountain  View,  CA  94043 

Society  of  Naval  Architects  and 
Marine  Engineers (Tech  Library) 

One  World  Trade  Center,  Suite  1369 
New  York,  New  York  10048 

Sun  Shipbuilding  &  Dry  Dock  Co 
Atin:  Chief  Naval  Architect 
Chester,  Pa  19000 

Sperry  Systems  Management  Division 
Sperry  Rand  Corporation (Library) 
Great  Neck,  N.Y.  11020 

Stanford  Research  Institute 
Attn:  Library 
Menlo  Park,  CA  94025 

Southwest  Research  Institute 
P.0.  Drawer  28510 
San  Antonio,  Texas  78284 
Attn :  App I i ed  Meehan i cs  Rev i ew 
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ane 

Texa^7872I 

Mr.  Roter^faggart 
94 1 1  Le®Cawy ,  Suite  P 
Fairfax,  V&V2203I 


^  «.  '  ■  ».  ».  ■.  4»,  T.  ^  Try.  »»  ■-  -y 


Ocean  Engr  Department 
Woods  Hole  Oceanographic  Inc. 

Woods  Hole,  Mass  02543 

Worcester  Polytechnic  Inst. 

A I  den  Reserach  Lab  (Tech  Library) 
Worcester,  MA  01609 

Applied  Physics  Laboratory 
University  of  Washington  (Tech  Library) 
1013  N.  E.  40th  Street 
Seattle,  WA  98105 

4  University  of  California 

Naval  Architecture  Department 
Berkeley,  CA  94720 
Attn:  Prof.  W.  Webster 
Prof .  J .  Pau I  I i ng 
Prof .  J .  Wehausen 
Library 

f  California  Institute  of  Technology 
Pasadena,  CA  91109 


Library 


Engineering  Research  Center 
Reading  Room 

Colorado  State  University 

Footh ills  Campus 

Fort  Collins,  Colorado  80521 


Florida  Atlantic  University 
Ocean  Engineering  Department 
Boca  Raton,  Florida  33432 
Attn:  Technical  Library 


Page  3 


GordonHteKay >Tbrary 
Harvard  uJHifersity 
Pierce  H?nl 
Cambridge,  MA  02138 

Department  of  Ocean  Engineering 
University  of  Hawaii  (Library) 
2565  The  Mali 
Honolulu,  Hawaii  96822 


Institute 
The  Uni 
Iowa  City, 
Attn:  Libr; 
Dr. 


raul  ic 
of  Iowa 
52240 


Research 


Prof.  Of  Phi Ufios 
Mechanics  MpaHmient 
The  John  Hopkins  UfH varsity 
Baltimore,  Md  21218 


Kansas  State  University 
Engineering  Experiment  Station 
Seaton  Ha  1 1 

Manhattan,  Kansas  66502 
Attn:  Prof.  D.  Nesmith 


St.  Anthony  Falls  Hydraulic  Laboratory 
University  of  Mi nnesota 
Mississippi  River  at  3rd  Avenue  S.  E. 
Minneapolis,  Minnesota  55414 
Attn:  m  1  "W^Ul 
Library 


/ 


Department  of  Naval  Architecture 
and  Marine  Engineering  -  North  Campus 
University  of  Michigan 
Ann  Arbor,  Michigan  48109 
Attn:  Library 


Davidson  Laboratory 
Stevens  Institute  of  Technology 
7 1 1  Hudson  Street 
Hoboken,  New  Jersey  07030 
Attn:  Library 


Applied  Research  Laboratory  Library 
University  of  Texas 
P.  0.  Box  8029 
Austin,  Texas  78712 


University  of  Kansas 

Chm  Civil  Engr  Department  Library  I 

Lawrence,  Kansas  60644 

Fritz  Engr  Laboratory  Library 
Department  of  Civil  Engr 
Lehigh  University 
Bethlehem,  Pa  18015 

Department  of  Ocean  Engineering 
Massachusetts  Institute  of  Technology 
Cambr i dge ,  MA  02139 
Attn:  Prof.  P.  Leehey 


Prof.  J.  Kerwin 


Engineering  Technical  Reports 
Room  1 0-500 

Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 


Stanford  University 
Stanford,  CA  94305 
Attn:  Engineering  Library 


Webb  Institute  of  Naval  Architecture 

Attn:  Library 

Crescent  Beach  Road 

Glen  Cove,  L. I .,  N.Y.  I  1542 

National  Science  Foundation 
Engineering  Division  Library 
1800  G  Street  N.  W. 

Washington,  D.  C.  20550 


Mr.  John  L.  Hess 
4338  Vista  Street 
Long  Beach,  CA  90803 

Dr.  Tuncer  Cebeci 
Mechanical  Engineering  Dept. 
California  State  University 
Long  Beach,  CA  90840 

Science  Applications,  Inc. 
134  Holiday  Court,  Suite  318 
Annapolis,  MD  21401 
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